Preview onto version 0.3.4

Here you can do some experiments with Rascal. The computation time is restricted to about 10 secs to prevent high load on server. Have a look at the documentation or at the examples. If you see wrong plots reload may help. You job is cancelled as soon someone else accesses this page, so consider a local installation.

Rascal-toolbox-0.3.4 (C) 2001,2002 Sebastian Ritterbusch

Enter 'help' for help or 'licence' for information about licence and warranty.

 "rascal.rc loaded, predefined e, i, pi"
>" Verification of the solution of linear (point) systems ";
>
>A=[1e+7 1;1e+7 1.0001];
>b=[2;7];
>R=1/A
 [   0.0010001000   -0.0010000000;-1.0000000000E+004 1.0000000000E+004]
>
>x=R*b
 [  -0.0049998000;5.0000000000E+004]
>A*x-b
 [7.2759576142E-012;7.2759576142E-012]
>
>" Iteration xk=(1-R*A)*xk+R*b starting with inflated x0:";
>xk=x*[0.8,1.2];
>i=0;
>C=(1-R*[1,1]*A);
>z=R*[1,1]*b;
>"[# Enclosure xk]";
>xn=C*xk+z;[i++ less(xn(1),xk(1))&less(xn(2),xk(2)) xk=xn]
 [1 1 [[  -0.0049998001,  -0.0049997999];[4.9999999999E+004,5.0000000001E+004]]]
>xn=C*xk+z;[i++ less(xn(1),xk(1))&less(xn(2),xk(2)) xk=xn]
 [2 1 [[  -0.0049998001,  -0.0049997999];[4.9999999999E+004,5.0000000001E+004]]]
>xn=C*xk+z;[i++ less(xn(1),xk(1))&less(xn(2),xk(2)) xk=xn]
 [3 0 [[  -0.0049998001,  -0.0049997999];[4.9999999999E+004,5.0000000001E+004]]]
>xn=C*xk+z;[i++ less(xn(1),xk(1))&less(xn(2),xk(2)) xk=xn]
 [4 0 [[  -0.0049998001,  -0.0049997999];[4.9999999999E+004,5.0000000001E+004]]]
>
>yk=(z=R*(b-A*[1,1]*x))*[0.8,1.2];
>C=(1-R*[1,1]*A);
>
>" Krawczyk Iteration yk=(1-R*A)*yk+R*(b-A*x) starting with inflated y0:";
>yk
 [[-1.7463171389E-014,2.6193447412E-014];[-2.6193447412E-007,1.7462298275E-007]]
>i=0;
>"[# Enclosure yk+x]";
>yn=C*yk+z;[i++ less(yn(1),yk(1))&less(yn(2),yk(2)) (yk=yn)+x]
 [1 1 [[  -0.0049998001,  -0.0049997999];[4.9999999999E+004,5.0000000001E+004]]]
>yn=C*yk+z;[i++ less(yn(1),yk(1))&less(yn(2),yk(2)) (yk=yn)+x]
 [2 1 [[  -0.0049998001,  -0.0049997999];[4.9999999999E+004,5.0000000001E+004]]]
>yn=C*yk+z;[i++ less(yn(1),yk(1))&less(yn(2),yk(2)) (yk=yn)+x]
 [3 0 [[  -0.0049998001,  -0.0049997999];[4.9999999999E+004,5.0000000001E+004]]]
>yn=C*yk+z;[i++ less(yn(1),yk(1))&less(yn(2),yk(2)) (yk=yn)+x]
 [4 0 [[  -0.0049998001,  -0.0049997999];[4.9999999999E+004,5.0000000001E+004]]]
>

Questions, Problems, Bugs? Please use this form.